function [Ae be] = expand(A, b)
be=[];
Ae = A;
[m n] = size(A);
for i = 1 : m
    be(i) = b(i) + 40*norm(A(i, :));
end
be=be';
end